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We perform the study of the stability of the Lorenz system by using the Jacobi stability analysis, 
or the Kosambi-Cartan-Chern (KCC) theory. The Lorenz model plays an important role for under¬ 
standing hydrodynamic instabilities and the nature of the turbulence, also representing a non-trivial 
testing object for studying non-linear effects. The KCC theory represents a powerful mathematical 
method for the analysis of dynamical systems. In this approach we describe the evolution of the 
Lorenz system in geometric terms, by considering it as a geodesic in a Finsler space. By associating 
a non-linear connection and a Berwald type connection, five geometrical invariants are obtained, 
with the second invariant giving the Jacobi stability of the system. The Jacobi (in)stability is a 
natural generalization of the (in)stability of the geodesic flow on a differentiable manifold endowed 
with a metric (Riemannian or Finslerian) to the non-metric setting. In order to apply the KCC the¬ 
ory we reformulate the Lorenz system as a set of two second order non-linear differential equations. 
The geometric invariants associated to this system (nonlinear and Berwald connections), and the 
deviation curvature tensor, as well as its eigenvalues, are explicitly obtained. The Jacobi stability 
of the equilibrium points of the Lorenz system is studied, and the condition of the stability of the 
equilibrium points is obtained. Finally, we consider the time evolution of the components of the 
deviation vector near the equilibrium points. 


* t.harko@ucl.ac.uk 
t chor.yin.ho@polyu.edu.hk 
t chun-sing-hkpu.leung@polyu.edu.hk 
§ stan.yip@polyu.edu.hk 



2 


I. INTRODUCTION 


Continuously time evolving dynamical systems are one of the basic theoretical tools for modeling the evolution of 
natural phenomena in every branch of physics, chemistry, or biology. Their usefulness in scientific/engineering applica¬ 
tions is determined by their predictive power, which, in turn, strongly depends on the stability of their solutions. Since 
in the measured initial conditions in a physical system some uncertainty inevitably does exist, a physically meaningful 
mathematical model must offer an understanding of the possible evolution of the deviations of the trajectories of the 
studied dynamical system from a given reference trajectory. Note that a local understanding of the stability is as 
important as the global evolution and control of late-time deviations. From a mathematical point of view the global 
stability of the solutions of the dynamical systems is described by the well studied theory of Lyapounov stability. In 
this approach the fundamental quantities are the Lyapunov exponents, measuring exponential deviations from the 
given trajectory It is usually very difficult to analytically determine the Lyapounov exponents, and therefore 

various numerical methods for their calculation have been proposed, and are used in various situations [^-[TTj|. On 
the other hand the local stability of solutions of dynamical systems is much less understood. 

Even that the methods of the Lyapounov stability analysis are well established, it would be interesting to study the 
stability of the dynamical system from different points of view, and to compare the results with the corresponding 
Lyapunov exponents analysis. Such an alternative approach to the study of the dynamical systems is represented by 
the so-called g:eometrodynamical approach, which was initiated in the pioneering work of Kosambi [l^ , Cartan [l3 | 
and Chern The Kosambi-Cartan-Chern (KCC) approach is inspired by the geometry of the Finsler spaces. Its 
basic idea is to consider that there is a one to one correspondence between a second order dynamical system and the 
geodesic equations in an associated Finsler space (for a recent review of the KCC theory see [H). The KCC theory 
is a differential geometric theory of the variational equations for the deviations of the whole trajectory to nearby ones 
M- In this geometrical description of the dynamical systems one associates a non-linear connection, and a Berwald 
type connection to the differential system, and five geometrical invariants are obtained. The second invariant, also 
called the curvature deviation tensor, gives the Jacobi stability of the system Il5l4l§l . The KCC theory has been 
applied for the study of different physical, biochemical or technical systems (see IITH^ L 

An alternative geometrization method for dynamical systems was proposed in and [23 | , and further investigated 
in [ 2 ^-[ 2 ^. Specihc applications for the Henon-Heiles system and Bianchi type IX cosmological models were also 
considered. In particular, in a theoretical framework devoted to a geometrical description of the behavior of 

dynamical systems and their chaotic properties was developed. 

In the Riemmannian geometric approach to dynamical systems one starts with the well-known results that the flow 
associated with a time dependent Hamiltonian H = 5°‘^PaPb/‘^ + V (a;“) can be reformulated as a geodesic flow in a 
curved, but conformally flat, manifold [l^. By introducing a metric of the form ds^ = W {x°') Sabdx°'dx^ , with the 
conformal factor given by W (x“) = E — V (x“), where E is the conserved energy associated with the time-independent 
iJ, it follows that the geodesic equation for motion in the metric gab = {x°') 5ab is completely equivalent to the 

Hamilton equations dx’^/dt = dH/dpa, dpa/dt = —dH/dxa [l^- This implies that the confluence or divergence of 
nearby trajectories a;“(s) and [a;-l-C]“('S) is determined by the Jacobi equation, i.e., the equation of geodesic deviation, 
which takes the form 


D2^a 

Ds^ 






( 1 ) 


where is the Riemann tensor associated with gab, and D/Ds = u°‘S/a denotes a directional derivative along 
= dx°'/ds. Linear stability for the trajectory x°'{s) is thus related to or, more exactly, to the curvature iF“. 
If, for example, everywhere negative, so that always has one or more negative eigenvalues, the trajectory 

must be linearly unstable |24l |. 

The global and local stability of solutions of continuously evolving dynamical systems was reconsidered, from a 
geometric perspective, in [30l| . It was shown that an unambiguous definition of stability generally requires the choice 
of additional geometric structures that are not intrinsic to the dynamical system itself. 

One of the most studied non-linear differential equations system is the Lorenz system which has been 

proposed as a possible description of the mechanism of the transition to weak turbulence in natural convection, by 
simulating thermal convection in the atmosphere, and modeling turbulent convective flow. An alternate view for the 
emergence of chaos in Lorenz-like systems was considered in |34l |. The effects of noise on the Lorenz equations in the 
parameter regime admitting two stable fixed point solutions and a strange attractor was studied in 1351 . Absolute 
periods and symbolic names were assigned to stable and unstable periodic orbits of the Lorenz system in |36| . Whether 
unstable periodic orbits associated with a strange attractor may predict the occurrence of the robust sharp peaks in 
histograms of some experimental chaotic time series for the Lorenz equations was investigated in [ 3 ^ . A nonlinear 
feedback approach for controlling the Lorenz equation was proposed in [s^. A particular return map for a class of 
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low dimensional chaotic models called Kolmogorov - Lorenz systems where studied from the viewpoint of energy and 
Casimir balance in [s^, by using a general Hamiltonian description. 

The possibility that the Lorenz system can be discussed in terms of the KCC-theory in Finsler space was pointed 
out in 1^, but without presenting any concrete analysis. A geometric viewpoint of the Lorenz system, based on a 
theory of tangent bundle, was proposed in [4l|. By introducing the geometrical viewpoints of second order system 
governed by Euler - Poincare equation or Lie - Poisson equation, the geometrical invariants of the Lorenz system have 
been obtained. It was shown that a torsion tensor, as one of geometrical invariants, relates to the chaotic behavior, 
characterized by the Rayleigh number, and results from the decomposition from the second order system to the 
tangent space (state space) and base space (configuration space), respectively. 

It is the purpose of the present paper to consider a full analysis of the Lorenz equations in the framework of the 
KCC theory. As a first step in this approach, the Lorenz system is reformulated as a set of two second order non-linear 
differential equations. The geometric invariants associated to this system (nonlinear and Berwald connections), and 
the components of the deviation curvature tensor are explicitly obtained, as well as its eigenvalues. The Jacobi stability 
of the equilibrium points of the Lorenz system is studied, and the condition of the stability of the equilibrium points 
is obtained. Finally, we consider the time evolution of the components of the deviation vector near the equilibrium 
points. 

The present paper is organized as follows. In Section |TT] we reformulate the Lorenz system as a set of two second 
order differential equations. The basics of the KCC theory to be used in the sequel are presented in Section IBII 
The Jacobi stability of the Lorenz system is considered in Section HVl The time evolution of the components of the 
deviation vector is studied in Section |Vl We discuss and conclude our results in Section IVIl 


II. THE LORENZ EQUATIONS 


In order to 
mechanics, we 


discuss the application of the KCC theory to systems of differential equations connected to fluid 
consider the Lorenz system of three nonlinear ordinary differential equations, given by 


I dX 
(j dt 


-X + Y, 


( 2 ) 


— = -XZ + pX-Y, (3) 

^ = XY-PZ, (4) 

where a, p and /3 are some free parameters. From a mathematical point of view these ordinary differential equations 
represent an approximation to a system of partial differential equations describing finite amplitude convection in a 
fluid layer heated from below. The Lorenz system results if the unknown functions in the partial differential equations 
are expanded in Fourier series, and all the resulting Fourier coefficients are set equal to zero except three [3l|. The 
parameters cr, p, and P can be interpreted from a physical point of view as the Prandl number, the Rayleigh number 
(suitably normalized), and the wave length number, respectively. Although the Lorenz set of equations is deterministic, 
its solution shows chaotic behavior for p > pcrit = o'((t + j3 + Y)/ {o — j3 — 1), and a > j3 + 1 [3l|, 


A. Second order differential equations formulations of the Lorenz system 

From Eq. m we can express Y as 


1 . 

Y = X + -X. 
a 


By substituting Y into Eqs. we obtain 

X + {1 + a)X + aXZ + cr (1 - p) A = 0. 
By taking the derivative of Eq. (U) with respect to the time we find 


1 •. 


Z = XY + XY-PZ = X[XY-X]+X[X + -X]- 


1 


X [X + ^x\ - pz 


( 5 ) 

( 6 ) 


( 7 ) 
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By substituting X as given by Eq. m we obtain the following equation for Z, 

Z+( + - - 2 ] XX - -X^ + (1-p +f3)X'^ + X'^Z - P'^Z = 0. 

\ a a J a 

At this moment we change the notation so that 

X = X^,X = Y^,Z = = Y^, 

and 

Y = X^, 

respectively. Hence the Lorenz system is equivalent with the following system of second order equations, 

J2 vl 

+ (1 + (t) + aX^X^ + cr (1 - p) = 0, 

+ - - 2 ) Alyl _ 1 ^ ^ ^ ^ ^^1^2 ^2 _ ^2^2 ^ Q 

dt^ \ (7 a J a ^ ' V / \ / 

Once and X‘^ are known, the variable X^ can be obtained as 

A^ = Ai + i-yi. 


( 8 ) 

(9) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


III. KOSAMBI-CARTAN-CHERN (KCC) THEORY AND JACOBI STABILITY 


In the present Section we summarize the basic concepts and results of the KCC theory (for a detailed presentation 
see M and [H). 

Let A4 be a real, smooth n-dimensional manifold and let TAl be its tangent bundle. On an open connected subset 
H of the Euclidian (2n +1) dimensional space i?" x i?" x i?i we introduce a 2n + 1 dimensional coordinates system 
(a;®,?/*,!), i = l,2,...,n, where (a;®) = (a;i,a;^,...,a;"), (?/®) = (?/i,...,y") and t is the usual time coordinate. The 
coordinates y® are defined as 


/ da;i dx'^ dx^'\ 


(14) 


In the following we assume that the time coordinate t is an absolute invariant. Therefore the only admissible 
coordinate transformations are 


t = t, i® = i® (a;i, a;^,..., a;") , i G {1, 2,..., n} . (15) 

Following [soil , we define a deterministic dynamical systems as a set of formal rules describing the evolution of points 
in a set S with respect to an external, discrete, or continuous time parameter t € T. More precisely, a dynamical 
system is a map [30[ 


cj) ■. T X S ^ S, {t, x) 1 -^ x), (16) 

which satisfies the condition (^{t, ■)o(j)(^s, •) = (j){t+s, •), Vt, s G T. For realistic dynamical systems additional structures 
must be added to this definition. 

In many situations the equations of motion of a dynamical system can be derived from a Lagrangian L via the 
Euler-Lagrange equations. 


d dL 
dt cly® 


dL 

9x® 


= F,,i = l,2,...,n, 


(17) 
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where i = 1,2,n, is the external force. The triplet {M, L, Fi) is called a Finslerian mechanical system [43l l44| . 
For a regular Lagrangian L, the Euler-Lagrange equations introduced in Eq. (El) are equivalent to a system of 
second-order ordinary (usually nonlinear) differential equations 


(Px^ 

dP 


+ 2G" {x\yPt) = 0,i G {1,2,...,n}, 


(18) 


where each function G* {xPyPt) is C°° in a neighborhood of some initial conditions ((x)q , (y)^ ,to) in fl. 

The basic idea of the KCC theory is to start from an arbitrary system of second-order differential equations of the 
form (1181) , with no a priori given Lagrangean function assumed, and study the behavior of its trajectories by analogy 
with the trajectories of the Euler-Lagrange system. 

To analyze the geometry associate to the dynamical system defined by Eqs. o, as a first step we introduce a 


nonlinear connection N on M, with coefficients iV), defined as [44 

N = - 

^ dyi 


(19) 


The nonlinear connection can be understood geometrically in terms of a dynamical covariant derivative 
for two vector fields v, w defined over a manifold M, we introduce the covariant derivative as 


V% = 


+ N{{x,y)w^ 


d 

dx' 


( 20 ) 


For Nf{x,y) = Tli(x)y^, Eq. (1201) reduces to the definition of the covariant derivative for the special case of a 
standard linear connection, as defined in Riemmannian geometry. 

For the non-singular coordinate transformations introduced through E qs. (1151) . we define the KCC-covariant differ¬ 
ential of a vector field {x) on the open subset fl C i?" x i?" x as [16l - [l9l | 


DC 

dt 


dt 


-N]C. 


( 21 ) 


For C = 2/* we obtain 


dt 


= Ny -2G^ = -e\ 


( 22 ) 


The contravariant vector field e* on il is called the first KCC invariant. 

We vary now the trajectories x*(t) of the system (IT8)) into nearby ones according to 

P (t) = x\t)+r]C(t), 


(23) 


where \ri\ is a small parameter, and ^*(t) are the components of a contravariant vector field defined along the path 
x'^jt). S ubstituting Eqs. (E^ into Eqs. (fTSl) and taking the limit y ^ 0 we obtain the deviation equations in the form 


[iMi 


Eq. 


d^C 

dC 




can be reformulate in the covariant form with the use of the KCC-covariant differential as 


dP 




where we have denoted 


5G* ; , 


dx'- 


/ 9m 

. WiV' H- - 

' ^ dt ’ 


(24) 


(25) 


(26) 


and we have introduced the Berwald connection G'^ji, defined as [1 

dN{ 

^ 

Qyl 


( 27 ) 
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Pj is called the second KCC-invariant or the deviation curvature tensor, while Eq. (I25p is called the Jacobi equation. 
When the system (IT^ describes the geodesic equations, Eq. (1^ is the Jacobi field equation, in either Riemann or 
Finsler geometry. . 

The trace P of the curvature deviation tensor is obtained as 


pini 

P = Pi = -2— - 2G'G* 

* dx^ 


y 


idNl 


-NlNl 


dm 


dx^ ' ' dt 

The third, fourth and fifth invariants of the system (1181) are defined as [if 


P® = 


1 


dP^ 

dy^ 


dy^ 


P® = ^ n® = ^ 

Qyl ^^Jkl- Qyl ■ 


(28) 


(29) 


The third invariant can be interpreted geometrically as a torsion tensor. The fourth and fifth invariants and 
^jki called the Riemann-Christoffel curvature tensor, and the Douglas tensor, respectively [H, [l^. In a Berwald 
space these tensors always exist. In the KCC theory they describe the geometrical properties and interpretation of a 
system of second-order differential equations. 

Alternatively, we can introduce another definition for the third and fourth KCC invariants, as 


where 


B;k = 


Sx^ 


_S_ 

(5x® 


d 

dx^ 


Sx^ ’ 

(30) 

- m—. 

® dy:i 

(31) 


The fourth KCC invariant can then be defined as 


O. _dBl, 

Bjki - 


dyi 


(32) 


In many physical, chemical or biological applications we are interested in the behavior of the trajectories of the 
dynamical system (I18|) in a vicinity of a point a;® (to). For simplicity in the following we take to = 0. We consider the 
trajectories a;® = a;®(t) as curves in the Euclidean space (P", (.,.)), where (.,.) is the canonical inner product of P®®. 
For the deviation vector ^ we assume that it obeys the initial conditions ^ (0) = O and ^ (0) = IT 7 ^ O, where O G P” 
is the null vector [l5l - [l8l| . 

Thus, we introduce the following description of the focusing tendency of the trajectories around to = 0: if ||C (^)ll < 
P, t ~ 0+, the trajectories are bunching together. On the other hand, if ||■C(^)|| > t^, t ~ 0+, the trajectories 
are dispersing [l5l - [l§ |. The focusing tendency of the trajectories can be also characterized in terms of the deviation 
curvature tensor in the following way: The trajectories of the system of equations (1181) are bunching together for t ^ 0'^ 
if and only if the real part of the eigenvalues of the deviation tensor PJ (0) are strictly negative. The trajectories are 
dispersing if and only if the real part of the eigenvalues of Pj ( 0 ) are strictly positive [iMl. 

Based on the above considerations we define the concept of the Jacobi stability for a dynamical system as follows 

[IM3: 

Definition: If the system of differential equations Eqs. m satisfies the initial conditions | |a;® (to) — i® (to)| | = 0, 
I [i® (to) — i® (to)| I 7 ^ 0, with respect to the norm ||.|| induced by a positive definite inner product, then the trajectories 
of Eqs. (flSl) are Jacobi stable if and only if the real parts of the eigenvalues of the deviation tensor Pj are strictly 
negative everywhere. Otherwise, the trajectories are Jacobi unstable. 

Graphically, the focussing behavior of the trajectories near the origin is shown in Fig. [1] 


The curvature deviation tensor can be written in a matrix form as 


with the eigenvalues given by 




A± 


1 

2 


Pi+P|±^ 


Pif + dP^ipf 


(33) 


(34) 
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FIG. 1. Behavior of the trajectories near zero. 


The eigenvalues of the curvature deviation tensor are the solutions of the quadratic equation 

A2 _ (pi + p|) A + - P^P^) = 0. (35) 

In order to obtain the signs of the eigenvalues of the curvature deviation tensor we use the Routh-Hurwitz criteria 
(45l| . according to which all of the roots of the polynomial P{X) are negatives or have negative real parts if the 
determinant of all Hurwitz matrices det Hj, j = 1,2, ..,n, are positive. For n = 2, corresponding to the case of 
Eq. (1551) . the Routh-Hurwitz criteria simplify to 

p^ + Pi < 0, Pi Pi - Pi Pi > 0. (36) 

A± describe the curvature properties along a given geodesic. Hence we can characterize the way the geodesic explore 
the Finsler manifold through the (half) of the Ricci curvature scalar along the flow, k, and the anisotropy 9, defined 

as [Hi 




p2 


and 


d = i(A+-A_) 


v/(pl_p2)2+4p^lp2 

2 


(37) 


(38) 


respectively. 


IV. JACOBI STABILITY OF THE LORENZ SYSTEM 

In the present Section we use the KCC approach for the study of the dynamical properties of the Lorenz system. We 
explicitly obtain the non-linear and Berwald connections, and the deviation curvature tensors for the Lorenz system. 
The eigenvalues of the deviation curvature tensor are also obtained, and we study their properties in the equilibrium 
points of the Lorenz system. The study of the sign of the eigenvalues allows us to formulate a basic theorem giving 
the Jacobi stability properties of the fixed points of the Lorenz system. 


A. The non-linear and Berwald connections, and the KCC invariants of the Lorenz system 


The Lorenz system can be formulated as a second order differential system, given by two equations of the form 


(PX^ 

dP 


+ 2G* {X\Y^) =0,i = l,2. 


(39) 









From Eqs. (El) and (fT^ it follows immediately that 


{X\X^, Y^) = 2 [(1 + 


and 


G^ (X\X^,Y^) = 


+ ^ - 2 ) x^Y^ - i (yi)" + (i-p + /3) (x^y + 


(x^fx^-/3^x^ 

respectively. Therefore we first obtain the components of the non-linear connection as 

, dGUX\X^,Y^) 1 ^ . 

- gyi =^(1 + ^),X,^=0, 


2 dG^(X\X^,Y^) 1/1 +a /3 

9yi ~ 2 [ a a 


N( = 


- 2 ] x^ —y\iv| = 0 . 

a 


For the components of the Berwald connection we obtain 


G 1 _ /^1 _ _ /^1 _ r \ /^2 _ /^2 _ /^2 _ /^2 _ f » 

11 — ^12 “ ^21 — ^22 — ^11 — ’ ^12 “ ^21 — ^22 ~ ^■ 

a 

The components of the first KCC invariant of the Lorenz system are given by 

+ aX^X^ + a(l - p)X\ 


and 


£2 = 


i + (i-p + /3) (x^f + (x^fx^ - /3^X^, 


respectively. 

The components of the curvature deviation tensor of the Lorenz system are given by 

pi = -crX^ -a(l-p) + ^(l + af, 

Pi = -aX\ 


P2 = 1 - 


2a 


yi 


1 — cr^ — 7/3(7 + (3 + 4(p — l)a 


4ct 


x^ - x^x"^. 


(40) 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


p^ = -{x^y + p\ 

For the trace of the curvature deviation tensor we obtain 

P = Pl+ Pi = - {Xy^ -aX^-a{l-p) + \{l + af + 

The time variation of the components of the deviation curvature tensor for the Lorenz system are represented in 
Figs. E1I51 

The third, fourth and fifth KCC invariants, as defined by Eqs. (1291) are identically equal to zero for the Lorenz 
system. However, the third invariant as defined by Eq. (1501) has a non-zero component. Generally, the third KCC 
invariant can be written as 


dN] dm dNi , dm, 

TDi _ k _|_ Arm ^ k _ ATI _^ 

~ dx>^ dX3 ^ dY^ dY^ ’ 


which has a single non-zero component. 


p 2 ^ridN? 1 + a 
Bi2 - N, . 


(52) 


(53) 


All the other KCC invariants of the Lorenz system are identically equal to zero. 
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FIG. 2. Time variation of the deviation curvature tensor component shown in the left figure, and the time variation of 

Piit), presented in the right figure, for cr = 10, p = 28, and /3 = 8/3, respectively. The initial conditions used for the numerical 
integration of the Lorenz system are X(0) = 1, T(0) = 5, and Z{0) = 10. 




FIG. 3. Time variation of the deviation curvature tensor component Pi{t), presented in the left figure, and the time variation 
of Piit), shown in the right figure, for a = 10, p — 28, and /3 = 8/3, respectively. The initial conditions used for the numerical 
integration of the Lorenz system are X(0) = 1, T(0) = 5, and Z(0) = 10. 


B. The Jacobi stability of the equilibrium points of the Lorenz system 


The three equilibrium points of the Lorenz system are S'o(0, 0,0), if p < 1, 

ip- 1)> {p- 1): P - 1 , and S'- 


— yj (3 {p — (3 {p — 1), p — 1 , if p > 1, respectively. From the point 

of view of the second order differential formulation of the Lorenz system and of the Jacobi analysis, the equilibrium 
points of the system given by Eqs. (ED and m are 


So{X^=0,X^ = 0),p<l, 


(54) 


and 


S4 


S_ 


Xl = ^(3{p-l),Xl = p-l\ ,p>l, 


X^_ = -Vl3ip-l),Xl=p-l 


,P>3, 


(55) 


(56) 


respectively. 

In the equilibrium points the components of the first KCC invariant vanish identically, so that 

e* (So) = P (S+) = e* (S_) = 0, i = 1,2. (57) 

We evaluate now the components of the curvature deviation tensor in the equilibrium points. First we obtain 


Pi (So) = -(1 - p)ct + 7 (1 + ct)\p < 1 


(58) 


PHSo)=P^{So) = 0,p<1 


(59) 
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For the equilibrium points 5”+ and S- we obtain 

Pi(5+) = i(l + a)2,p>l, 


and 


P 2 (5'+) = -cr\//3(p- 1), p > 1, 


52 /c ^ _ Vl^ip- 1) (- 7 / 3 cr + /? - 0-2 + 


P{iS+) = 


4a 


> 1 , 


P^{S+) = p^-/3ip-l),p>l, 


PliS_) = -il + af,p>l, 


P 2 {S-) = ay/j3{p- 1),p > 1, 


22 /a ^ _ V^iP- 1) (/5(7o' - 1) + cr^ - 1) 


(*5-) = 


4a 


p > 1 , 


P|(5_) = / 32 -/ 3 (p- 1 ),p> 1 , 


respectively. 

The eigenvalues of the curvature deviation tensor in the equilibrium points are obtained as 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


('S'o) — ^ [o’(4p + (T — 2) + 1], A_ (S'o) — /3^, p < 1, 


(69) 


A+ (5+) = 


12 


/3(/?-p+l)--(cr + l)2 


+ /3(p - 1) [/3(7o- - 1) + 0-2 - 1] + 


/3(/3 — p + l) + -(cr + 1)^ kp> 1, 


A_(5+) = -^ - 



/3 (/3 - p + 1) - -(cr + 1)2 + /3(p - 1) [/3(7cr - 1) + cr2 - 1] + 


/3 (/3 — p + 1)) + —-—^l,p > 1, 


(70) 


(71) 


A+ ^S.) = -< 


/3(/3 - p + 1) - -(cr + 1)2 


+ /3(p — 1) (/3(7cr — 1) + 0-2 — 1) + 


/3(/3 — p+1) + -(ct + 1)^ C)P> 1) 


A-(^-) = 2^- 



n 2 


/3(/3-P+1)--(ct + 1)2 


+ /3(p - 1) (/3(7cr - 1) + cr2 - 1) + 


(72) 


/3(/3 — p + 1) + - (cr + 1)^ >, p > 1. 


(73) 
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The eigenvalues of the deviation curvature tensor have the property 

A+ (5+) = A+ (5_), A_ (^+) = A_ (S.) . (74) 

The trace k of the deviation curvature tensor, as well as the anisotropy 0 of the Lorenz system are obtained as 

^ (‘S'o) = 2 4 [o'(4p + O' — 2) + 1]| , (75) 

^ (‘5'o) = 2 + ^^(dp + cr - 2) + 1]| , (76) 

/c {S+) = K {S-) = i [4/3(/? - p + 1) + (a + 1)^] , (77) 

and 


9{S+) = 0{S.) 


/3(/3-p+l)--(cr + l)2 


+ /3(p - 1) [/3(7u - 1) + cr2 - 1], 


(78) 


respectively. 

Taking into account the previous results we can formulate the following theorem, giving the Jacobi properties of 
the equilibrium points of the Lorenz system: 

Theorem, a) The equilibrium point Sq (0, 0) of the Lorenz system is Jacobi unstable. 

b) If the free parameters /3, p > 1, and cr of the Lorenz system satisfy simultaneously the constraints 

/3{P-p+l) + ^{a + lf <0, (79) 


and 


{/3 [-7pcr + p + cr(cr + 9)] - 2cr(p - l)(cr + 1)} > 0, 


(80) 


respectively, then the equilibrium points S'+ P ~ 1 

Jacobi stable, and Jacobi unstable otherwise. 


and S 4 . 


— \/l3[p — 1), p — 1 of the Lorenz system are 


V. THE ONSET OF CHAOS IN THE LORENZ SYSTEM 

The behavior of the deviation vector i = 1, 2, giving the behavior of the trajectories of a dynamical system near 
a fixed point (to) is described Eqs. (l24ll and (1251) . In the case of the Lorenz system these equations can be written 
generally as 

—+ (cr + 1)~J^ + O' [(1 - P) + ^ 2 ] = 0, (81) 

and 

+ [(/3 - o + 1) - 2Yi] + [2a(/3 - p + X 2 + l)Xi + (/3 - a + l)Yi] ^^(t) + 

a{X,-/3){l3 + X,)e{t)=0, (82) 

respectively. The deviation vector is obtained from its components as 


m = \/[e{t)f + {e{t)f 


(83) 
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In order to obtain a quantitative description of the onset of chaos in the Lorenz system, we introduce, in analogy 
with the Lyapounov exponent, the instability exponents 5i, z = 1, 2, and <5, defined as 


SAS) = lim - In 

' t->oo t 


^iO 


,* = 1 , 2 . 


and 


S{S) = lim - In 




t^oc t I ^10 

In the following we investigate the behavior of the solutions of Eqs. dni) and 
Lorenz system. 


(84) 

(85) 

near the critical points of the 


A. Behavior of the deviation vector near 50(0,0) 

Near the equilibrium point 50(0,0) the deviation equations Eqs. (1811) and (15^ take the form 


dt^ 


dt 


and 


dt'^ 


-/3T(<)=0, 


( 86 ) 


(87) 


respectively. In this case the deviation equation for the origin of the Lorenz system can be separated in two independent 
equations, with the general solutions given by 


e{t) 

ClO 


and 


i(v^4po7j:(o71T)2_o--l)t ^ gi(--\/4pCT+(cr-l)2-o'-l)i 

\/4pcr + (cr - 1)2 

o-0t (e2/3t _ 1) 


e{t) = 


2/3 


?20, 


( 88 ) 


(89) 


where we have used the initial conditions ^^(0) = 0,^^(0) = ,^io, and ^^(0) = 0,^^(0) = ^ 20 , respectively. The time 
behavior of ^^(t) is determined only by the coefficient /3 of the Lorenz system. For the deviation vector we obtain 

■jlosinh-O) 


m = 


10 


4/32 


4pcr + (cr — 1)2 


1/2 


(90) 


The time dependence of the deviation vectors and is represented, for different values of the parameters cr, p, 
and /3, in Fig. |31 

For a fixed /3 and cr, with the increase of the parameter p, the deviation curvature component increases very 
rapidly in time, indicating the onset of chaos in the Lorenz system. In the large time limit the increase of ^^(t) is 
exponential. The deviation curvature vector component ^2 also increases exponentially in time, a behavior that is 
independent of the values of /3. The time variation of the absolute value of the deviation vector ^(t) is represented in 
Fig. [S] 

With the increase of p the absolute value of the deviation vector rapidly increases in time, indicating the onset of 
chaos in the Lorenz system. With the help of Eqs. (1551) and (15^ we immediately obtain 


(5o) = 


1 


4pcr + (cr — 1) — cr — 1 


and 


<52 (5o) = /3, 

respectively. The instability exponent 5 can be estimated as 


1 ! p 

5 1 In < ^20 e 


2Pt 


o-\/4pcr+(cr-l)2t 


2t \ 4/32 4pCT + (cr - 1)2 


(91) 

(92) 

(93) 


The time variation of the instability exponent 6 is represented, for fixed values of the parameters P and cr, in Fig. [51 



























13 




/: 


/ 

/ ■ 

/ / 
/ 


/ 

/ 

/ / 


/ / 


0.0 0.1 


I I I I I I I I I I I I I I I I I I I I I I I I I I 1 I 

0.2 0.3 0.4 0.5 0.6 0.7 



FIG. 4. Time variation of the deviation vector components (left figure) and ^^(t )/^20 (right figure), in the vicinity of 

the equilibrium point So, for different values of the parameters /3, p and a. In the left figure a = 10, /? = 8/3, and p = 15 (solid 
curve), p = 20 (dotted curve), p = 25 (dashed curve), and p = 28 (long dashed curve), respectively. In the right figure /3 = 2/3 
(solid curve), /3 = 5/3 (dotted curve), /3 = 8/3 (dashed curve), and P = 11/3 (long dashed curve), respectively. 



FIG. 5. Time variation of the absolute value of the deviation vector ^(t)/^io for /3 = 8/3, cr = 10, ^lo = 10~®, ^20 = 10~®, and 
for different values of p: p = 10 (solid curve), p = 20 (dotted curve), p = 25 (dashed curve), p = 28 (long dashed curve), and 
p = 33 (ultra-long dashed curve, respectively. 



FIG. 6. Time variation of the instability exponent 5 for P = 8/3, cr = 10, ^10 = 10~®, ^20 = 10~®, and for different values of p: 
p = 10 (solid curve), p = 20 (dotted curve), p = 25 (dashed curve), p — 28 (long dashed curve), and p = 33 (ultra-long dashed 
curve, respectively. 
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B. 


Dynamics of the deviation vector near S'+ ^j3(p — 1), p — 1 and S_ —\//3(p — 1), p — 1 


For both fixed points 5"+ y^/3(p — 1), p — 1 


and S- 


1))P- 1 


dynamics of the deviation vector near the given fixed points take the form 


the differential equations describing the 


^+ 1 )^+= 0 , 

and 

respectively. 

The behavior of the components of the deviation curvature vector is shown in Fig. [T] 


(94) 


(95) 




FIG. 7. Time variation of the deviation vector components C^(t) (left figure) and ^^{t) (right figure) in a logarithmic scale 
near the equilibrium points S± /3(p — 1)^, for P = 8/3, a = 10, and for different values of p: p = 15 (solid curve), p = 20 
(dotted curve), p = 25 (dashed curve), p = 28 (long dashed curve), and p = 33 (ultra-long dashed curve, respectively. The 
initial conditions used to integrate the deviation equations are ^^(0) = ^^(0) = 0, C^(0) = 10“^°, and 5^(0) = 10~^, respectively. 

The time variation of the instability exponent S (5+) = 6 (S-) is represented in Fig. [S] 



FIG. 8. Time variation of the instability exponent S (S+) = <5 (S_) near the equilibrium points S± ^±-\//3(p — 1) j, for P = 8/3, 
cr = 10, and for different values of p: p = 15 (solid curve), p = 20 (dotted curve), p = 25 (dashed curve), p = 28 (long dashed 
curve), and p = 33 (ultra-long dashed curve, respectively. The initial conditions used to integrate the deviation equations are 
^^(0) = ^^(0) = 0, 5^(0) = 10“^°, and C^(0) = 10“®, respectively. 
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C. The curvature of the deviation vector 


In order to obtain a quantitative description of the behavior of the curvature deviation tensor in the following we 
analyze the signed geometric curvature kq of the curve ^{t) = , which we define, according to the standard 

approach used in the differential geometry of plane curves as 


kq{S) = ko = 





+ 


iHt 


i,ji ■ 


(96) 


If we denote = a and cr + 1 = 6, we obtain an explicit formula for 


Ko, 


Ko = 


'?loC20 


oi-b/2)t 


8a 


+ (a + + 6“-®*]^! 

(a - b){2p -a + + (a - b){b -2/3- + 

(a + 6)(2/3 + a + + (a + b){a + b - 2^)e(-“/2-/5)‘ 


3/2 


(97) 


In the following we restrict our study to the equilibrium point iS'o(0,0), and we fix the values of the parameters as 
a = 10, /3 = 8/3, ^10 = 10“^°, and ^20 = 10“®, respectively. The variation of the curvature kq is represented in Fig. [HI 



FIG. 9. Time variation of the curvature kq of the deviation vector ^{t) near the equilibrium points So (0,0), for j3 = 8/3, 
CT = 10, ^^(0) = 10~^°, and ^^(0) = 10“®, respectively, and for different values of p: p — 10 (solid curve), p = 15 (dotted curve), 
p = 20 (dashed curve), p = 25 (long dashed curve), p — 28 (ultra-long dashed curve, and p = 33 (ultra-ultra-long dashed 
curve), respectively. 


As one can see from the figure, for the chosen range of the parameters the curvature of the two-dimensional deviation 
vector is positive for small values of time, reaches the value zero at a certain moment Iq, enters the region of the 
negative values, and then it tends to zero in the limit of large times. Near the time origin 1 = 0 the curvature of the 
deviation vector curve can be represented in a form of a power series as 


Ko{t,p) 


11 4545vToTp-|-45584vToTi , ll(877185\/ToTp-l-2748853\/ToT) 2 , 

202\/ToT 9272709 ^ 1873087218 ^ 


(98) 


In the hrst approximation the time interval to for which kq (to, p) = 0 is given by 


9999 

2(4545p-h 45584) 


1.099 

p-h 10.02' 


(99) 


Therefore the time interval after which the curvature kq changes sign is an indicator of the development of chaos 
in the Lorenz system. There exist therefore a critical time interval so that if the transition from positive to 
negative kq occurs at times t > tg’'**, the evolution of the Lorenz system is predictable, and deterministic, while if 
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the transition occurs for time intervals t < then the underlying evolution of the Lorenz system will be chaotic 
in the long term. The curvature becomes zero when the two components of the deviation vector satisfy the condition 

{to) = e (to). 


D. Comparing with linear stability analysis 

We will now compare Jacobi stability and linear stability at the equilibrium point S'o(0,0). The linearized Lorenz 
system at (0,0, 0) is given by 


X' 


—a a 

0 


'a' 

Y 

= 

p -I 

0 


Y 

Z 


0 0 

-/3 


Z 


We see that the equation for Z is decoupled and Z{t) tends to 0 exponentially fast. The time dynamics of X and 
Y is governed by the two-dimensional system 


' X' 


—a 

a 


a' 

Y 


P 

-1 


Y 


Let A := 


—a fj 

P -1 


. Let T = trace(A) and A = det(A). Then, 


T + \/t‘^ — 4 A r — — 4A 

-^1 =-- 1 A2 = — 


2 2 
On the other hand, in Jacobi stability analysis, we find the eigenvalues of Tj (0, 0) are 

A-h (-S'o) = i [cr(4p + a-2) + l] = - 4A), 

A_ (5o)=/3^P<l, 


( 100 ) 

( 101 ) 


Thus, one of eigenvalues of PJ(0,0) recovers some information of the stability of the system at the origin. That is, 
the origin is a center or spiral, when A+ (S'o) < 0. 


VI. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the stability analysis of the Lorenz system from the point of view of the 
KCC theory, in which the dynamical stability properties of dynamical systems are inferred from the study of the 
geometric properties of the Finsler space geodesic equations, equivalent with the given system. By transforming the 
Lorenz to an equivalent system of two second order differential equations, the dynamical system can be interpreted 
as representing the geodesic motion of a ’’particle” in an associated Finsler space. This geometrization of the Lorenz 
system opens the possibility of applying the standard methods of differential geometry for the study of its properties. 
We have obtained, and analyzed in detail, the main geometrical objects that can be associated to the Lorenz system, 
namely, the non-linear connection, the Berwald connection, and the first, second and third KCC invariants. The main 
result of the present paper is the Jacobi stability condition of the equilibrium points of the Lorenz system, showing 
that the origin is always Jacobi unstable, while the Jacobi stability of the other two equilibrium points depends on 
the values of the parameters of the system. By considering the standard values of cr, (3 and p in the Lorenz system, 
a = 10, P = 8/3, and p = 28, it turns out that for this choice of the parameters all the equilibrium points are Jacobi 
unstable. From the point of view of linear stability analysis, for 0 < p < 1, the zero fixed point S'o = (0, 0,) is globally 
stable; from a physical point of view this refers to the non convective state. For p > I and p not too large, the state 
with roll convection, referring to the S± equilibrium points, is stable From a physical point of view this means 
that the phase space consists of two regions, which are separated by the stable manifold of the zero fixed point; the 
trajectories which start in one of the two regions are attracted by the corresponding nonzero fixed point. 

We have also considered in detail the behavior of the deviation vector near the equilibrium points. In order to 
describe the behavior of the trajectories near the equilibrium points we have introduced the instability exponent 6, 
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as well as the curvature kq of the deviation vector trajectories. The curvature kq of the curve ^(t) can be related 
directly to the chaotic behavior of the trajectories via its transition moment from positive to negative values. An 
early transition indicates the presence of chaotic states. Therefore we suggest the use of the curvature of the deviation 
vector as an indicator of the onset of chaos in non-linear dynamical systems. In it was suggested that the torsion 
tensor geometrically expresses the chaotic behavior of dynamical systems, i.e. a trajectory of dynamical systems with 
the torsion tensor is not closed. By using the same definition as in it turns out that indeed there is a non-zero 
torsion tensor component = —(1 + cr)/2(7 ^ 0,\la ^ 0. Therefore the existence of chaos in the Lorenz system is 
intimately related to a non-zero a, which is needed for the chaotic behavior of the Lorenz system [s^. 
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